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Abstract 



Molecular dynamics results for the dynamic structure of a Heo.77Neo.23 gas 
mixture at two densities (15.8 nm~'^ and 36.1 nm~'^) are presented. A clear 
description of the crossover from hydrodynamic modes to distinct excitations 
for the two species is obtained. The higher density data neatly show high- 
and low-frequency branches setting on in the dispersion curve, with a rather 
localized transition. The lower density results remarkably agree with exist- 
ing neutron scattering data and, differently from previous simulation studies, 
display hydrodynamic behavior up to A; ~ 2 nm~^. A smooth transition to 
fast sound is shown to take place for 2 < k/nroT^ < 5, where the present MD 
data fill the existing gap in the experimental results. 



PACS numbers: 61.20.Ja, 61.20.Lc, 78.35.+C 



1 



Typeset using REVT^ 



Inelastic neutron scattering studies of the dynamic behavior of two-component fluids 
were usefully performed on dense He-Ne and He-Ar gas mixtures Both conventional 

and, recently, Brillouin p neutron spectroscopy produced, in particular, a rather 
extended set of data for the dynamic structure factor S{k, u) of He-Ne mixtures with low 
concentration (20-30%, typically) of the heavier component, and total number density n 
around 16 nm~^. Important and detailed indications on the nature of dynamic excitations 
in these mixtures have been obtained from the above neutron data. On the one hand, they 
show that at wavevectors k > 4 nm~^ the total S{k,uj) contains two Brillouin doublets, 
one at frequencies ±uji^\ corresponding to a "fast" sound propagating with velocity ci nearly 
equal to that of pure He at the same temperature and total density of the mixture (~ 540 
m s~^), and one, not giving rise to visible inelastic features in S{k,u!), at lower frequencies, 
±ujs , attributed to the presence of a "slow" sound related to the Ne atoms dynamics (as 
in previous papers, labels 1 and 2 refer to He and Ne, respectively). On the other hand, 
from high resolution neutron Brillouin spectra at k < 1.8 nm^^ |^ no decoupling of the 
dynamics of the two species is observed. Instead, one single acoustic mode was detected, 
with a linear dispersion in remarkable agreement with the hydrodynamic sound velocity 
c ~ 370 m s~^ of the mixture, as estimated from the van der Waals equation 0. There is, 
therefore, experimental evidence that the transition from hydrodynamic behavior towards 
distinct excitations for the two components occurs in the wavevector region around /c ~ 3 
nm~^, for this He-Ne state. 

This experimental result, however, does not agree with the analysis of a molecular dy- 
namics (MD) simulation of the same mixture 0, where the onset of a lower frequency mode 
in the Ne dynamics was detected at a fc- value as small as 0.6 nm~^ by analyzing the partial 
correlations. At larger k values, the MD data of Ref. indicate a progressive departure 
from hydrodynamic sound (more evident in the Ne case), characterized by split, low- and 
high-frequency, branches for the heavy and light component, respectively. Very recently, 
new simulation results on the same mixture were also announced which, while agreeing 
with the experimental data of Ref. [^j for the total S{k, u), again confirm the findings of the 
older simulation. 

In order to investigate this apparent discrepancy between neutron Brillouin and MD 
results, we performed extended simulations of two Heo.77Neo.23 mixtures: one in the same 
thermodynamic state [n = 15.8 nm~^, T = 39 K) as in the quoted neutron and MD 
I^J^ studies; the other in a far denser state {n = 36.1 nm~^, T = 60 K), because, with 
increasing density, hydrodynamics is expected to hold in a wider k range, and the features 
of the Rayleigh-Brillouin spectrum should be enhanced by the reduced line widths and by 
the higher speed of sound in the mixture, both effects allowing for a clearer detection of 
dynamic decoupling of the two components. 

The constant- A^V^i? simulations were carried out by using 87808 particles confined in a 
cubic box of 17.7 nm side at the lower density (13.4 nm side at the higher density) with 
periodic boundary conditions. The equations of motion were integrated using the leap-frog 
algorithm with a time step of 20 fs at low density and of 12.5 fs at high density. After 10^ 
time steps for thermalization, the dynamics has been followed for about 2400 ps and 600 ps, 
at low and high density, respectively. The time step for configuration recording was fixed 
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to 0.1 ps. The corresponding wavevector ranges were 0.4 < fc/nm~^ < 5, at low density, 
and 0.5 < k/nm~^ < 6.6, at high density. The He-He and Ne-Ne interactions were modeled 
according to the Aziz et al. potentials plpl, while the He-Ne cross interaction parameters 



were derived by applying the Lorentz-Berthelot rules ||10 . 

Here we first summarize the high density simulation results. At low k, i.e. in the 
hydrodynamic regime, each Sij{k,uj) (z, j = 1, 2) is predicted to be composed of a central 
line, resulting from the superposition of two Lorentzian components (both located at zero 
frequency) related to concentration diffusion and heat diffusion, and of two Lorentzian lines, 
with asymmetric correction, centred at ± a;s, i.e. the typical Brillouin doublet PJTI|. In 
formulae: 



S{k, uj) 



TT 



OOI^OI Q02^02 Qs^s + &s(^ + ^s) fls^s - bs{u - UJs 



Other models based on extended sets of dynamical variables have been applied in the- 
oretical studies of He-Ne mixtures However, the eight-parameter fit of Eq. (1), 



convoluted with the simulation resolution function, to the MD data provides an excellent 
description of all Sij{k,uj) in the whole /c-range of our simulations. In the upper frames of 
Fig. 1, Sii{k,u!) and S22{k,u!) are plotted at three k values, together with the fitted lines. 
Logarithmic scales are used to enhance the visibility of the inelastic features. The frequency 
range in the plots is limited to three times the fitted ooi^^ , though the used fit range is larger. 
We note that Eq. (1) proves very good over more than 4 orders of magnitude in intensity. 
The lower frames show both data and fitted curves multiplied by u;^, i.e. the longitudinal 
current spectra Cij{k,uj) = u'^Sij{k,u), including the cross term Ci2{k,u). The following 
features can be observed. The Brillouin peak in Su{k,uj) is visible at all k, while that of 
S22{k,uj) appears only as a very weak shoulder at the intermediate k, and can no longer be 
seen at the highest one. Cij{k,u) presents another structure at a lower frequency, of the 
order of the width of the central line of Sij{k, u), i.e. where interdiffusion processes dominate 
the spectrum. The main feature of Cii{k,uj) is, at all k, the Brillouin peak, while in the 
Ne-Ne correlation the intensity evolution of the two peaks as a function of k determines the 
position of the maximum of C22{k,u): at low A;, this is close to ujf'\ but at higher fc, where 
the Brillouin intensity becomes too low to produce a maximum, it shifts to the location of 
the lower frequency peak. Inspection of the cross correlation reveals that, in the first case, 
the maximum of C22{k,uj) occurs where Ci2{k,Lj) is positive, indicating in-phase motions 
of the two components, while, in the second case, it is located in a frequency range that 
cannot be related to an acoustic propagation, since it is characterized by negative values of 
Ci2{k,uj), a signature of inter-species diffusion dynamics. 

The dispersion curves ujs{k) of the fitted inelastic peak frequency in Sii{k,u!) and 
S22{k,ui) are shown in Fig. 2. Up to A; ~ 4 nm~^, the He and Ne curves clearly show 
the same (substantially linear) behavior with a slight positive dispersion: the collective ex- 
citations propagate with speed Ci = C2 = c varying from 900 to 1000 m s~^. Above this k 
value, the dispersion curve clearly splits into two branches, indicating the onset of a dynamic 
decoupling of the two species that gives rise to fast and slow sound in the mixture, with 
ci = uj^^ /k rising up to 1050 m s~^ and C2 = ujf'' /k decreasing sharply down to 800 m s~^, 
at = 6.5 nm~^. These values are in reasonable agreement with the speed of sound of the 
pure components (at the same temperature and total density of the mixture), i.e. 1150 m 
s-i for He [|l2l and 840 m s'^ for Ne 
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Fig. 2 also shows the dispersion curve derived as the position of the maximum in 
Cij{k,uj), an often used procedure In the He case such a determination provides a 

frequency larger than uji^\ while for the Ne-Ne correlation this method produces, for the 
reason explained before, a wrong dispersion curve which departs from linearity at rather 
small k values, with frequencies much lower than uos , and in the characteristic range of the 
concentration diffusion processes. The wiggles in such a curve, evident in Fig. 2, correspond 
to a range where the two peaks in C22{k, uj) have comparable heights, so that its maximum 
jumps from one position to the other due to statistical noise. 

We now turn to the low density simulation results, for which we repeated the partial 
structure factors analysis described above, obtaining, on the whole, a similar picture. Again, 
Eq. (1) remarkably accounts for all Sij{k, u) at all k values, as shown in Fig. 3. The inelastic 
features of each Cij{k, u), including the sign of the cross term, display a k dependence similar 
to that of the high density mixture. However, in the Ne case the Brillouin peak amplitude 
vanishes for k > 2 nm~^ and prevents from a reliable determination of its position. Thus, 
the Ne-Ne inelastic structures cannot be followed further. 

The fitted uj^^^ and cjp'* as functions of k are given in Fig. 4. No departure from 
hydrodynamic sound is observed for either component (ci = C2 = c ~ 370 m s~^) for all k 
values below 1.9 nm~^. For 2 < k/nm~^ < 5, the propagation speed Ci of the light particles 
excitations neatly undergoes a smooth transition from the hydrodynamic sound speed c of 
the mixture to the quoted fast sound velocity (~ 540 m s~^). Since the total S{k,u!) is 
dominated by the Su{k,u!) term and shows therefore inelastic peaks at frequencies similar 
to those of the He-He correlation, this result provides the missing link between the low- and 
high- A; neutron data (see Fig. 2 of Ref. ^ and Fig. 3 of ref. 

The dispersion curves obtained from the maxima of the longitudinal current spectra, 
also shown in Fig. 4, again give a different indication: they both would suggest a departure 
from the hydrodynamic straight line at k around 0.7 nm~^, with a deviation reaching 15% 
around k ^ 2 nm~^ for He [l^ and much larger (negative) values for Ne. Indeed, for the Ne 
case, such a curve is in good agreement with the analogous one in Fig. 4 of Ref. 0. 

The molecular-dynamics investigation reported in this Letter clearly illustrates the onset 
of decoupled modes in dense gaseous binary mixtures, from which the fast sound phenomenon 
originates. The He-Ne mixtures here considered provide valuable insight on these effects. In 
the high density case, the splitting of the dispersion curve into two branches is a sharp effect 
and a slow sound propagation is clearly detected. At the lower density, a very good agreement 
is found with the existing neutron measurements. The present MD data smoothly bridge 
the gap between the low-fc neutron Brillouin results ^ and the older high-/c experiment [|l|. 
In particular, hydrodynamic behavior up to A; ~ 2 nm~^ is confirmed, and a continuous 
transition to fast sound is shown to take place for 2 < A;/nm~^ < 5. From the present data 
it cannot be decided whether a lower branch of the dispersion curve also exists at the lower 
density, similarly to what happens at the higher one. It may be that the Ne dispersion curve 
follows that of He up to a /c value larger than 2 nm~^, before the splitting takes place, or 
that no acoustic mode propagates through the Ne atoms above k ^ 2 nm~^, as suggested 
by the vanishing amplitude of the Brillouin peak. Finally, it is clear from this work that, at 
least in the presence of very weak inelastic structures compared with quasi-elastic diffusion 
processes, an analysis based on the location of maxima in the longitudinal current spectra 
can lead to quite incorrect results. 
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FIGURES 



FIG. 1. MD spectra at the higher density and three values of k. Upper frames: Sii{k,uj) (dots) 
and S22{k,uj) multipHed by 3 (circles). Lower frames: Cii(A;,a;) (dots), C22ik,u;) multiplied by 3 
(circles), and Ci2{k,uj) (squares). In all frames the solid lines are best fits with Eq. (1). Only few 
data points are shown for the sake of clarity. Error bars, where not visible, are smaller than the 

(2) 

size of the symbols. The arrows in the lower frames show the values of . 

FIG. 2. Dispersion curves of the acoustic modes in MD Sn(k,uo) and (A;, w) at the higher 
density. The fitted Wg^^ (dots) and Wg^^ (circles) are shown with error bars. Solid lines are obtained 
as positions of maxima of Cii(A;,a;) (upper) and C22ik,Lo) (lower). 

FIG. 3. MD spectra at the lower density and three values of k. The format of the picture is the 

(2) 

same as in Fig. 1. The arrow in the bottom right frame is missing since cjg cannot be obtained 
reliably from the fit for k > 2 nm~-^. 

FIG. 4. Dispersion curves of the acoustic modes in MD Sii{k,Lj) (top) and S22{k,to) (bottom) 
at the lower density. The fitted Wg^^ (dots) and Wg^'' (circles) are shown with error bars. Solid lines 
are obtained as positions of maxima of Cii(A;,a;) (top) and C22(A;,a;) (bottom). Dashed straight 
lines correspond to propagation speeds of 374 and 544 m s~^. 
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